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If three right-handed neutrinos are added to the Standard Model, then, for the three known 
generations, there are six quarks and six leptons. It is then natural to assume that the symme- 
try considerations that have been applied to the quark matrices are also valid for the lepton 
mass matrices. Under this assumption, the solar and atmospheric neutrino data can be used 
to determine the individual neutrino masses. Using the x 2 fit, it is found that the mass of the 
^ ■ lightest neutrino is (2-5) x 1CT 3 eV, that of the next heavier neutrino is (10-13) x 1CT 3 eV, 

<D \ while the mass of the heaviest neutrino is (52-54) x 10~ 3 eV. 
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I. INTRODUCTION 

In the work of Lehmann, Newton and Wu []lp, the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix is expressed in terms of the masses of the three generations of quarks: 



d) is) lb 



(1.1) 



This is accomplished by introducing a new horizontal symmetry. Some of the earlier attempts in 
this direction are given in [|3], U, while some of the more recent work on this topic is to be found 
in [|5], ||], for example. 

Recent experiments at Super-Kamiokande [|7], [8[] indicate the presence of neutrino oscillations, 
which would imply that the neutrinos are not all massless. If it is accepted that the neutrinos are 
not massless, then it is most natural in the Standard Model to introduce three right-handed 
neutrinos in addition to the three known left-handed ones. In this way, there are six quarks and 
six leptons. For recent reviews of neutrino physics, see JT0|] . Most of the recent work on neutrino 
masses is focused on the ideas of grand unification and the see-saw mechanism (see, e.g. [pi]]). 

In this paper, the consequences of a universal quark-lepton mixing are studied. While many 
authors favor Majorana masses for the neutrinos, this universality is most natural if they instead 
have Dirac masses. 1 In other words, the method of [JTp is used to express the lepton CKM matrix 
and the neutrino mixing matrix in terms of the masses of the three generations of leptons: 

v P \ fv„\ (v., (12) 



el \ /i j \t 

Of course the masses of the three charged leptons are accurately known, leaving as unknown 
parameters the masses of the three neutrinos. Thus there are three parameters to be determined 
instead of seven, the three masses plus the four in the lepton CKM matrix. 

It is the purpose of this paper to use the data from solar neutrinos [T3], [R], [15], [T7p 
and atmospheric neutrinos []7p to determine the three neutrino masses separately, not only the 
differences of their squares. As compared with earlier work on this model [|T%{], mixing in the 
charged lepton sector is also taken into account. 

In Sees. H and[II| we review the case of quarks and apply the model to the case of leptons. In 
Sec. |V|, the rotation matrix is discussed, and in Sec. [V] we construct the resulting mixing matrix for 
charged-current interactions. In Sec. |VT|we review the relevant formulas for atmospheric neutrino 
propagation, whereas Sees. |V1I| and |V111| are devoted to the three-flavor MSW problem and solar 
neutrinos. In Sec. |IX] we collect the results of the fits to the data, and Sec. |X| contains a discussion. 



Some properties of the mass matrix are discussed in Appendix |A|, whereas Appendices |B] and |C] 
contain technical details of the analytical solution of the MSW equations. 



1 As shown in [|l2p, neutrino-oscillation experiments cannot distinguish between massive Majorana and Dirac neu- 
trinos. 
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IL REVIEW FOR THE CASE OF QUARKS 
For the quark mass matrices, the result of Ref. [jlj] is 

d(d) \ 

M(d) = | d(d) c{d) b{d) , (2.1) 
b(d) a\d)j 

id(u) \ 
M(u) = | -id(u) c{u) b{u) , (2.2) 
b(u) a(u)J 

with 

b 2 (d) = 8c 2 (d), b 2 {u) = 8c 2 (u). (2.3) 

The diagonalization of these mass matrices is achieved by the orthogonal matrices R(d) and 
R(u), explicitly 2 

M(d) = R(d)diag(m d ,-m s ,m b )R T (d), (2.4) 

M(u) = diag(i, 1, l)R(u)diag(m u , — m c , m t )R T (u)dmg{— i, 1,1). (2.5) 

In terms of R(d) and R(u), the CKM mixing matrix [0] is written as 

V = J R T (w)diag(-i, 1, l)R(d). (2.6) 

For both the u quarks and the d quarks, the number of independent parameters in the mass 
matrix is three. Hence they can be expressed in terms of the three quark masses. The relations are 

a + c = Si — m 3 — m 2 + mi, 
8c 2 + d 2 — etc = —S 2 = m 3 m 2 — m^mi + m 2 m 1 , 

ad 2 = —S 3 = m 1 m 2 m3, (2.7) 

where mi, m 2 , m 3 mean respectively m u , m c , m t and m^, m s , m& for the w and d quarks. Because 
of the known masses, the inequalities 

rrii — m 2 — m 3 (2.8) 



are used throughout the analysis of Ref. O. 
In Ref. [|Tp , it is found that [JT^] 



\J CP \ ~ 2.6 x 10" 5 . (2.9) 



This is consistent with the experimental value of (3.0 ± 1.3) x 10 5 [p0|]. This experimental value 
is expected to improve significantly in the near future. 



Note misprints in Eqs. (27) and (35) of Ref. [|l[j, 
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III. APPLICATION TO THE CASE OF LEPTONS 

Since the neutrinos are now known to have masses, the symmetry considerations for the quark 
mass matrices are equally applicable to the leptons. Furthermore, the leptons, not being confined, 
have masses that are better defined than those of the quarks. Replacing d and u by the charged 
leptons £ and the neutrinos v, Eqs. fl2.ip, ( |2.2D and Q2.3D take the form 



d{£) 

M{£) = | d(£) c{£) b{£) | , (3.1) 
b{£) a(£) j 

id(u) 

M(u) = | -id(u) c{v) b{v) | , (3.2) 
b{v) a\v) 

with 

b 2 (£) = 8c 2 {£), b 2 {v) = 8c». (3.3) 
The rotation matrices R{£) and R(v) are defined in exactly the same manner: 

M(£) = R(£)diag(m e ,-m fM ,m T )R T {£), (3.4) 

M{y) = diag(i, 1, l)i?(i/)diag(mi, -m 2 , m 3 )i? T (z/)diag(— i, 1, 1), (3.5) 

where in Eq. ( p.5[ ) the masses of the three neutrinos are designated as m l5 m 2 and m 3 . The lepton 
CKM mixing matrix is 

V 1 = i? T (z/)diag(-i, 1, 1)R{£). (3.6) 
There are actually significant differences between the quark case and the lepton case. 

(A) In the quark case, that there is CP violation has been known for many years [ [2T| |22| ]. In 
the lepton case, it is not known whether CP is conserved or not. While it is tempting, on 
the basis of quark-lepton universality, to believe that CP non-conservation also holds for 
leptons, the possibility of lepton CP conservation cannot be excluded. In the former case, 
Eq. ( p.2[ ) holds; in the latter case, both i and — i on the right-hand side of Eq. Q3.2| ) are 
replaced by 1. To take both possibilities into account, Eq. (CO) needs to be generalized to 



ed{y) 

M{y) = | e*d(v) c(v) b(v) \ , (3.7) 
b{v) a\u) 

where the two cases above correspond to e = i, 1, respectively. Similarly, Eqs. fl3.5D and 



( pT6[ ) take the forms 

MO) = diag(e, 1, l)R(v)&ag(m u -m 2 , m 3 )i?»diag( e *, 1, 1), (3.8) 

and 

V 1 = i? T (z/)diag(e*, 1, 1)R{£). (3.9) 
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(B) In the case of quarks, both the masses and the absolute values of the elements of the CKM 
mixing matrix are known experimentally. In the case of the leptons, there is much less 
experimental information. While the masses of the three charged leptons are accurately 
known, the corresponding knowledge about the neutrinos is limited to differences between 
the masses squared, mf , mf , m\. It is for this reason that the lepton case may be considered 
to be more challenging, and one of the first tasks is to determine the individual masses from 
the existing experimental data. This is to be carried out in Sec. 

(C) The lack of knowledge about the neutrino masses has another profound consequence. For 
quarks, the inequality ( [0] ) holds for both u and d. In contrast, because of the presence of 
the minus sign with m 2 in Eq. Q3.5D, all that is known about neutrino masses is 



m i < m 3- (3.10) 

In other words, for neutrinos, Eq. ( |3.10| ) can be used, but not (|2TE|). The first task is therefore 
to determine the allowed region in the space (mi,m 2 ,m 3 ) of neutrinos, which must be 
between those permitted by (|2.8[) and (Oil). 



It follows from ( [277D that the parameter a, which must be positive, satisfies the cubic equation 

9a 3 - 17S ia 2 + (8S? + S 2 )a - S 3 = 0. (3.11) 
Any real cubic equation can have either one or three real solutions. Where there is one real 



solution, that one is negative, and thus unphysical, as is seen from (|2.7D. Where there are three 
real solutions, one of them is negative, while two are positive. We shall refer to these two positive 
solutions as Solution 1 (larger a) and Solution 2 (smaller a). 

It is instructive to consider briefly the simple case m\ = m 2 = 0. In this case, it follows from 
(|Z77[) that Si = m 3 and S 2 = S 3 = 0, and ( |3"TTD reduces to 



9a 3 - 17m 3 a 2 + 8m 2 3 a = 0, (3.12) 

with the solutions 

a = m 3 , ~m 3 , 0. (3.13) 

Here is the limiting value of the negative solution and hence is of no interest. From the above 
definitions, in this case a = m 3 is Solution 1 while a = | m 3 is Solution 2. From (|3.7|), the mass 
matrices are 

M{u) = | ] m 3 , (3.14) 

for Solution 1, and 













m 3 , 


V 
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M{v) = i ^ m 3 , (3.15) 



for Solution 2. 

In Ref. [[U] for the quark mass matrices, only Solution 1 was considered. 

The above considerations on the cubic equation ( |3.11| ) can be used to determine the allowed 
physical region in the (mi/m 3 ,m2/m 3 ) plane, as shown in Fig. |T]. This region is only slightly 
larger than the triangle given by the inequality (pT8|), with two additional regions, one where m 2 > 
m 3 and the other a very small one with m\> m 2 - 

The following notation is convenient. 
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1.5 




FIG. 1 : The allowed region for the three neutrino masses mi , m-2 and m 3 is bounded by parts of the axes 
and the solid curves C\ and C 2 . The lines S\ and 52 correspond to mi = m-2 and m-2 = m 3 , respectively. 

(1) 1Z denotes the entire allowed region, not including the boundary, for the three neutrino 
masses m x , m 2 and m 3 , as shown in Fig. [Tj. 

(2) IZo denotes the part of 1Z where mi < m 2 < m 3 . 

(3) TZi denotes the part of 1Z where m 1 > m 2 . 

(4) IZ2 denotes the part of 1Z where m 2 > m 3 . 

(5) Si denotes the part of 1Z where m x = m 2 ; similarly S 2 the part where m 2 = m 3 . 
Thus 

n = n Q + n 1 + n 2 + s 1 + s 2 . (3.16) 

(6) C\ is the curved part of the boundary of 1Z\. 

(7) C 2 is the curved part of the boundary of 1Z 2 . 

(8) C3 is the plane 

< m 2 /m 3 < (17+ ~ 1.42, m x = 0. (3.17) 

(9) C 4 is the plane 

< m l /m 3 < 17- 12v^ - 0.029, m 2 = 0. (3.18) 

IV. THE ROTATION MATRIX 

There is an elementary but somewhat complicated issue of the sign ambiguities in the definition 
of the rotation matrices R{v) and R(£) as given by Eqs. (^^) and (|3~3|). 
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It is seen from Eq. ( |2.7D that, for either Solution 1 or Solution 2 as defined in Sec. the values 
of m 1 , m 2 and m 3 in 1Z determine those of a, b 2 and d 2 , but not b and d. Therefore, for either R{u) 
or R(£), there are actually eight distinct _R's in 1Z: R k±± , where k = 1 for Solution 1 and k = 2 
for Solution 2, and where the first ± and the second ± designate the "parities" (or sign factors) of 
b and d respectively (see below). 



Since these i?'s are determined by 



r 


d 






d 


c 


b 


\r = r\ 


\o 


b 


a) 






(4.1) 



the elements of R are given explicitly by 

R n = d(mx - a)yi, R 12 = -d(m 2 + a)y 2 , Ri 3 = d(m 3 - a)y 3 , 

R 21 = m 1 (m 1 - a)y 1 , R 22 = m 2 (m 2 + a)y 2 , R 23 = m 3 (m 3 - a)y 3 , 

R 31 = bmxyi, R 32 = -bm 2 y 2 , R 33 = bm 3 y 3 , (4.2) 

where the values of yi, y 2 and y 3 are such that R is orthogonal. 

This condition of orthogonality does not determine the signs of these y's: 

±1 

Vj = , =, (4.3) 
yd 2 (\j - a) 2 + A|(Aj - a) 2 + b 2 \ 2 

where 

Ai = mi, A 2 = -m 2 , A 3 = m 3 . (4.4) 



The entire problem is to choose the three ± signs in Eq. (|4.3|). 

Strictly speaking, any choice of sign will do. Since such choices lead to a large number of 
possible CKM mixing matrices for leptons (for each R, there are 2 3 = 8 possible choices of 
signs), it is useful to make an intelligent choice of these signs. 

The basic principle to be used to choose the signs is that of continuity, i.e., the continuity of the 
nine R^ for each R. For example, each of these i?^ must be continuous in 1Z. Note that the con- 
tinuity of an R implies the continuity of its derivatives with respect to the masses m 1 , m 2 and m 3 . 

The problem to be solved is to find a set of eight 3x3 matrices R k±± with the following 
conditions: 

(1) R k±± are continuous in 1Z; 

(2) tfjf* > in K ; and 

(3) R 1 ^ = R 2±± ond. 

Condition (3) requires the following comments. First, this condition can be imposed on C\ or 
C 2 , but not on C 3 or C 4 . The reason is that, on C\ and C 2 , the values of a for the first and second 
solutions are the same. Secondly, what this equation means is that, given an R 1 , for example 
R 1++ , there is a choice for the two ± signs for R 2±± so that the i?'s are equal; this can only be 
achieved on either C\ or C 2 , not on both. Actually, there is no such choice: If C 2 is chosen, no 
solution exists because of the condition (2). It is therefore necessary to impose this condition on 
Ci, for any "6 parity" and "d parity." This is already indicated in condition (3) above. 
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TABLE I: Signs of the coefficients of the cubic equation (4.5) 
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S? + S 2 


—Sl£»2 + ^3 


n 2 
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± 
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+ 




± 


+ 




+ 




+ 





TABLE II: Signs of c for Solution 1 and Solution 2. 





Solution 1 


Solution 2 


n 2 










+ 


Tit 


+ 


+ 



The "parities" or sign factors are defined as follows: For positive "6 parity," b has the opposite 
sign of c, which is the convention of [P. In particular, for Solution 1, b is then positive in TZo. 
When the "6 parity" is flipped, R^ 3 , R 23 , -R 3 i and R 32 change sign. The "d parity" is the sign of d 
in all of 7Z. When this sign is flipped, R\ 2 , R\ 3 , R 21 and R 31 change sign. 



Attention is next turned to the cubic equation Q3.12| ) for a. Since c = Si — a from Eq. ([2.7D, the 
corresponding cubic equation for c is 

9c 3 - 10S lC 2 + (S 2 , + S 2 )c - (SiSz - S 3 ) = 0. (4.5) 

The nice formula 

S3 - StS 2 = (m 2 - mi)(m 3 - m 2 )(m 3 + mi) (4.6) 

shows that S 3 — SiS 2 is positive in 7Z , and in fact implies that, in 7Z , two of the solutions c from 
Eq. ( |4.5p are positive, while one is negative. In a similar way, the signs of the coefficients of the 
cubic equation ( ft~5D are listed in Table |. 



Since the three solutions of this cubic equation are known to be all real in 7Z, Table follows 
immediately from Table |. In particular 

c = on (Si for Solution 1 ; 

c = on S 2 for Solution 2. (4.7) 

The next task is to show that 

a - mi > (4.8) 
in 7Z, where the equality sign holds only on S 2 , and then only for Solution 2. Similarly, 

m 3 - a > (4.9) 

in 7Z, where the equality sign holds only on <Si, and then only for Solution 1. These derivations 
are straightforward and hence omitted. Note that the left-hand sides of these inequalities play an 



important role in Eq. (fO)). The two relations may conveniently be summarized as 

mi < a < m 3 . (4.10) 

cS 2 ,Sol2 .Si.Soll 
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Solution 1 + + 




FIG. 2: Solution 1: (is) matrix elements vs. m\ and m,2 (wi3 = 1), for z = 1,2,3 and j = 1, 2, 3. The 
dark regions indicate where the values exceed 0.5, and the lightly shaded regions indicate where the values 
are between and 0.5. In the unshaded (white) regions, the values are negative. 



It is also useful to note that 

iSi,Soll £2, Sol 2 

but in this latter case, these are not extrema. 

With this knowledge, the resulting signs are easily determined as given in Figs. ^ and |3|. The 
signs of yj, respecting conditions (1), (2) and (3) above, is such that: 

• For Solution 1, some Rjj vanish on <Si. They are ij = 13, 23, 31 and 32, corresponding to 
m 3 — a and b becoming zero. (These thus have opposite signs in 1Z and 1Zi.) 

• For Solution 2, some -R^ ±=t vanish on c>2. They are ij = 11, 21, 32 and 33, corresponding to 
a — mi and b becoming zero. (These -Rj 2 ? ±=t thus have opposite signs in 1Z and TZ 2 -) 
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Solution 2 + + 




FIG. 3: Solution 2: i?? ++ (V) matrix elements vs. m\ and ra 2 (mg = 1). Compare with Fig. ^| 



In order to have Rfj that are continuous, it is required to flip signs of some yj as these boundaries 
Si and S 2 are crossed. 

The mixing matrices display strong variations with mi and m 2 . For Solution 1, the diagonal 
elements dominate in much of the parameter space, whereas for Solution 2, this is not the case. 

For the charged leptons, the masses are strongly hierarchical. Thus, the rotation matrices cor- 
respond to the lower left-hand corners of those displayed in Figs. ^| and |[ For Solution 1, this is 
close to the unit matrix, whereas for Solution 2 certain non-diagonal elements are also significant. 
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V. THE MIXING MATRIX 

We write the unitary (but not necessarily real) mixing matrix as 







u e2 


u e3 


u = 












u r2 


u r3 



where (cf. Eq. 

U = (V e y = R T (e) diag(e, 1, l)R{u) (5.2) 
relates the neutrino mass eigenstates to the flavor states: 

\u e ) = U e i\vi) + U e2 \u 2 ) + U e3 \u 3 ), (5.3) 

etc. These are the states which enter in charged-current interactions. 

We shall label these mixing matrices by the Solutions (1 or 2) involved in the rotation ma- 
trices R(v) and R(£), and the "parities" associated with the parameters b and d, as discussed in 
Sec. |V]. Actually, only the product of the "^-parities" and that of the "d-parities" matter. Thus it 
is convenient to define 

b parity = (b parity)^ x (& parity)^, 

d parity = (d parity)^ x (d parity)^. (5.4) 

Therefore, there are a total of 

2x2x2x2x2 (5.5) 

different U matrices, where two factors of 2 arise from the two Solutions for the v and t sectors, 
two factors of 2 arise from the "6 parities" and the u d parities," and the last factor of 2 arises from 
the choice of CP conservation or CP non-conservation (e = 1 or e = i). 

For a representative case, Solution 2 for R(u) and Solution 1 for R(£), b parity — —,d parity = 
+, and no CP violation, we show in Fig. |] the resulting mixing matrix [7(2, 1) — h Since, for the 
charged leptons, Solution 1 is close to the unit matrix, the resulting U(2, 1) — h is rather similar 
to the corresponding R{v) for Solution 2 and parities — h 



VI. ATMOSPHERIC NEUTRINOS 



Let us next review the oscillation formulas that are relevant for atmospheric neutrinos. A 
neutrino state, which was a pure v a state at t = 0, will at time t have the form 



«(*) = U <*i e iEjt ^ ( a = e ' V' r ) 

3=1 



(6.1) 



For relativistic neutrinos, the energy is given as Ej ~ p + m3/2E, and the survival probability for 



the takes the form 



P„ 



1 -4 



IM 2 |CU 2 sin 2 



+ \U&\ \U^\ z sin 



2 „:„2 f ^ m 32^ 



l^irl^srsm 2 ' 



AE 



V AE 



(6.2) 
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FIG. 4: Mixing matrix [U (2, 1) — h]ij vs. mi and 7712 (TO3 = 1), for i = 1, 2, 3 and j = 1, 2, 3. The dark 
regions indicate where the values exceed 0.5, and the lightly shaded regions indicate where the values are 
between and 0.5. In the unshaded (white) regions, the values are negative. 



where Amf- = mf — m?. 

After propagation over some distance, a neutrino of a different flavor may appear. The proba- 
bility amplitude for the transition a — > (3 is given by 

(Mt)Wa(0)) = Y, U k U ^e l < t/2E . (6.3) 

k 

In contrast to the survival probability, this expression is not invariant under complex conjugation 
off/. 
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A case of particular interest is 

= \ Y, u ** u ^ eimlt/2E \ 



(Am 2 f\ 1 
-^f-J | + cyclic, (6.4) 



where we have used the orthogonality, J2k U* k U[ik = 0. For the real case, this simplifies to 



4E 



^UrsU^UnU^ sin 2 ( — ^ ) . (6.5) 



It is instructive to study the simple limit 

|A7t4l| < l Am 32l> kmlyt/AE < 1, (6.6) 
with all [/ a & = 0(1). Then, by unitarity, Eq. ( |63| ) simplifies to the familiar expression 



P„(t) ~ 4|^ 3 | 2 |f/.3| 2 sin 2 {^-) ■ (6.7) 

Fitting the data within a two-flavor model, with sin 2 9 and Am 2 as independent parameters, 
one finds [0] large mixing angles. In the limit of Eq. (|6~7o|), this corresponds to large values for 
|tV7r3| 2 ,seeEq.f77D. 

The observed suppression of atmospheric [0] suggests masses of the order m ~ 0.05 eV. 
In order to determine the neutrino masses, we formulate a x 2 by comparing predicted z/ M and v e 
fluxes with data: 

. erf 

The experimental data used are those from Super- Kamiokande [[7J]: the 8 data points (bins in 
E/L, where L = ct) for and the 8 data points for u e . These sixteen data points are treated as 
separate inputs, allowing an overall normalization constant for the two sets of data. Also, since 
the various survival and transition probabilities are rather sensitive to the precise values of energy 
and oscillation length, we averaged over these, within each of the 8 bins. 

We show in Fig. |] the contributions to x 2 from the atmospheric-neutrino data, for the mixing 
matrices U (2, 1) — + corresponding to CP conservation (left part) and CP non-conservation (right 
part). The figure shows x 2 a s a function of mi and m2, for fixed m 3 = 0.05 eV. 

The different solutions and parities that determine the mixing matrices give rather differently 
shaped x 2 minima when plotted vs. m\ and m 2 . For most cases, the minima occur inside the region 
TZq. For others, they occur near Si or near <S 2 - 

Comparing with Figs. |2| and 0, we see that Solution 1 provides large mixing for m 2 being a 
sizable fraction of m 3 , whereas Solution 2 favors relatively smaller values of m 2 , or m 2 close to 
mi. 
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VII. THE THREE-FAMILY MSW MECHANISM 



The coupled equations satisfied by the three neutrino wave functions are []2 

'M 2 i Mf 2 M 2 3 
Ml M| 2 M| 3 
Ml Ml Mf 3 , 









(Mr)\ 


)] 


1 


02 (r) 






\Mr)J 



(7.1) 



where D(r) = y/2G-pN e (r), with G F the Fermi weak- interaction constant and N e (r) the solar 
electron density at a distance r from the center of the sun. Furthermore, we denote the flavor states 
v e = <f>i, = </> 2 , v T — 3 . These are the states which enter in charged-current interactions. 
The evolution of the neutrino wave functions is determined by the squared mass matrix, 




ecd 
b 2 + c 2 + d 2 
b(a + c) 




'Mil 
Ml 

Ml 



M 2 2 
M 2 2 
Ml 



M? 3 
Ml 
ML 



(7.2) 



the neutrino momentum, p, and the solar electron density. Here, M 2 - = (M 2 )y , and e = 1 (CP 
conservation) or e = i (CP non-conservation). 

It is actually a good approximation to take an exponential electron density, N e (r) = 
N e (0) exp(— r/r ). A fit to the solar density as given by [ |27| ] gives r = 6.983 x 10 4 km. For this 
case of an exponential solar density, the three-component wave equation can be solved in terms of 
generalized hypergeometric functions, 2F2 p4]]. 



The case treated in [f24p was that of a real mass matrix. In that case, by scaling and shifting the 
radial variable, u = r/r + uq, with u determined such that 



D(0)r e uo = 1, 
Eq. dTTTl ) could be transformed into the form 

. d 
du 



(7.3) 



>1 


»1 




'ux + e~ u 


X2 


^2 






X2 




>3 






. X3 








^3. 





>1 






^2 


:«) 




.^3 


;«)J 



(7.4) 



with cui, cj 2 , u> 3 , X2 and %3 all real. 

We now have to address a small complication due to the possible non-reality of the mass matrix 
induced by CP non-conservation, and the fact that also the charged lepton states have to be rotated. 
Consider the case of CP non-conservation, i.e., e = i. In order to follow as closely as possible the 
procedure of [123]], we need to rotate to a neutrino basis which in the absence of matter (D(r) = 0) 
becomes that of the mass eigenstates. This now involves diagonalizing the lower right-hand part 
not of [MO)] 2 , Eq. 0, but of 



N 2 = R T (£)M 2 R(£) 



(N 



2\ 

11 

2^\* 

12 
2^* 
13 



(iV 2 ) 12 
(N 2 



22 



(N T 



113 



23 



(N 2 )h (iV 2 )33. 



(7.5) 



where not only the (1,2) and (1,3) elements are complex [together with (2,1) and (3,1)] (cf. 
Eq. Q7.2D), but also the (2,3) and (3,2) elements. This prevents a simple diagonalization like in 
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However, we may make the lower right-hand part real by a unitary 2x2 rotation U. In other 
words, we perform the transformation 



U 



\N 2 ) 22 (N 2 ) 23 





x y 


, with 


fa 


= u 


fa 




y z_ 




.fa. 




.fa. 



(7.6) 



By also defining 



cos 6*o 


— sin 6q 


r 


X 


V 




sin 9 


cos 6> 


2p 


JJ 


z 





together with 



r ° (N 2 ) 



n, 



and 



2p 

fa(u) = fa(u) 



X2 


_ ?~0 


.X3_ 


~ 2p 



.^3(w) 



cos 6*o sin 9q 
— sin cos 9 



cos #o — sin 6 

sin# cos^o 

cos 6*o — sin # 

sin #0 cos 6*o 





UJ 2 


o" 







o> 3 _ 



(7.7) 



(iV 2 ' 

L(iv 2 ; 



12 



'13. 





fa{u) 




.fa («)_ 



(7.8) 



(7.9) 



then Eq. ( f7TT| ) can be written in the desired form (J7J3]). The a/s are real, but the \2 and \3 will in 
general be complex. The rotations among fa and 3 given by U and # need not concern us, since 
we are here only interested in the electron neutrino, v e . 

When w 2 ^ W3, these fa, fa and fa can be expressed uniquely in terms of a single fa 



fa 
fa 
fa 



. d 

% U2 

du 



d 



X2 I * 

X3 (i 



d 

du 
d 

du 



U3J 
^2 J ip, 



(7.10) 



where tp satisfies the third-order ordinary differential equation [ |24| ] 



d 



du 



l- /il U M2 



d 



du 
d 

l- LU 2 

du 



(7.11) 



Here, //i, /z 2 , and /i 3 are the eigenvalues of the right-hand matrix in Eq. ( |7.4| ) (without the term 
e~ u ), ordered such that 

fJ>l < V>2 < fJ>3- (7.12) 

Equation ( p.l 1| ) is the differential equation for the generalized hypergeometric function 2 F 2 — 
see, for example, p. 184 of p5|]. Three linearly independent solutions of this third-order differen- 
tial equation ( T7 . 1 1| ) are 
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if, 



(i) 



K l e~ i ^ u 2 F 2 



1 - z(^ 2 - A*l' 



1 - i(/i 3 - /ii) 



(2) 



(3) 



K 2 e~^ u 2 F 2 



K 3 e-^ U 2 F 2 



1 - z(/ii - /i 2 ), 1 - «(a*3 - A*2y 



-2(^2 -A*3j, 
1 - z(^i - /i 3 ) 



~H^3 - yU3j 
1 - i(ji2 - fi 3 ) 



le 



le 



le 



(7.13) 



where K, K\, K 2 and K 3 are arbitrary non-zero constants. Since Eq. < J7 . 1 ID is linear, the general 
solution is 

i> = C^ (1) + C 2 ^ + 0^ (3) , (7.14) 

from which the ipi, ip 2 and ^3 can be obtained using Eq. ( |7.10| ). 

For the case of two flavors, the products in Q7.1 1| ) consist of one less term each, and a famil- 
iar confluent hypergeometric function 1 Fi (also denoted Whittaker function or parabolic cylinder 



function) is obtained [fZq]. 

These functions are trivial when u — > oo. In fact, outside the sun, they can be approximated by 
the exponential pref actors, since 



2-^2 



0*1, 



a 2 
b 2 







1. 



(7.15) 



In order to impose the boundary conditions that only electron neutrinos are produced in the 
sun, we have to determine these functions at large and negative values of u. The series expansion 
is in principle convergent, but it is not practical for large absolute values of both parameters and 
the argument. One possible way of dealing with these generalized hypergeometric functions has 



been given in []Z4]]. The procedure used there is as follows. First, ip\ , ip\ and ip\ J ' are evaluated 



,(3) 



approximately using Barnes' integral representation for 2 F 2 [p5[|. Since we have not managed to 
apply this same procedure to the other ifj's, they are expressed in terms of another generalized 
hypergeometric function 3 F\. Since 3 Fi has an integral representation in terms of the usual hyper- 
geometric function 2 F\, these 3 Fi can be evaluated by numerical integration. The choice of the 
contours of integration has been discussed in detail in p4]]. 

For completeness we give in Appendix [B] asymptotic formulas for these 3 Fi. These asymptotic 
formulas turn out to be quite useful and in particular are accurate for the region of the minimum 



2 to be discussed below. 



X 

Some details on the book-keeping of reconstructing the neutrino wave functions from the 2 F 2 
and 3 Fi are given in Appendix |C]. 



VIII. SOLAR NEUTRINOS 



In order to compare the predictions of the model to data, we form a x 2 by comparing the predic- 
tions to the available flux data. For the solar-neutrino flux, we take the values given by the 'BPOO' 
solar model p7|]. For the solar-neutrino data, we use the total rates from the Chlorine experiment 
[]T3[], the Gallium experiments [ |14| , p3| ] (we average the two results), the Super- Kamiokande ex- 



periment [|Sp , and the SNO experiment [|T^|, |T7]]. We adopt the neutrino energy spectra and detector 
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efficiencies as given by Bahcall et al. [|28|], and, for Super- Kamiokande and SNO, we also include 
the neutral-current cross section O]. We do not consider the day-night effect, since this is con- 
sistent with zero [jSp. Neither do we consider the electron-recoil spectrum, since this is consistent 
with being flat p^| ] (see, however [^Op). For the solar flux, we integrate over the spectrum (E u ), 
taking into account the detector efficiency e(E v ): 

$ = I dE » ^A E M E v)Pu{E v ). (8.1) 
i J 

In Fig. |5| we show the contributions to x 2 from the solar-neutrino data, as functions of m x and 
m 2 , for m 3 = 0.05 eV. As opposed to the atmospheric-neutrino data, the solar-neutrino data give 
a minimum x 2 that is well localized in the mi-m 2 plane, with little dependence on m 3 . 

The lower panels in Fig. [5| give the corresponding total x 2 = xltm + Xsoiar Since the minimum 
in xLar 1S ram er well-localized ("steep"), the location of the overall minimum (for fixed m 3 ) is 
largely determined by the minimum in Xsoiar- However, the extent to which the minimum in xltm 
overlaps with that of x 2 olar determines how good the fit is. 

The cases considered here are {7(2, 1) — +, with CP conservation (left) and CP non- 
conservation (right). Similar scans for the other mixing matrices (other Solutions, other parities) 
yield rather similar results. There is little difference between the cases of CP conservation and 
CP non-conservation. 



IX. RESULTS 
A. Combining Atmospheric and Solar Data 



Starting at the minima found from two-dimensional scans as described in Sees. |VT] and |VIII| , 
but now with mi, m 2 and m 3 all free, we find total minima for a given solution of the mixing 
matrix. The fitting has been performed using two different procedures. Procedure A is the method 
proposed by Hata and Langacker [jnj, where we allow for correlations between different data. 
Procedure B is a more transparent approach, where the data are treated as uncorrelated, and the 
SNO data are left out. There are 16 degrees of freedom. The detailed numerical results differ 
somewhat between these two procedures. 

The best such results from fitting procedure A are collected in Table |D. Two cases stand out. 

They are Solution "(2, 1) " and "(2, 1) - +," with x 2 values ranging from 16.3 to 18.1. Thus, 

Solution 2 is favored for the neutrino mass matrix, whereas Solution 1 is favored for the charged 
lepton mass matrix. Also, negative "6-parity" is favored, whereas there is no clear preference for 
a particular "d-parity." The latter observation is commented on in Sec. 

With fitting procedure B, these two Solutions have x 2 ranging from 15.8 to 17.0 as can be seen 
from Table |TV|. 

There is a strong clustering of m 3 values at 0.052-0.054 eV, for different Solutions and with 
different "parities." However, there are also a few fits which are not much inferior with m 3 values 
around 0.04 eV, and one at 0.095 eV. These have a x 2 that is higher by about 4 units. 



The best-fit mass values may roughly be related to the conventional Q32| , p3Q atmospheric and 
solar neutrino parameters as Am^ ~ ml — m\ = 2.6 x 10 -3 eV 2 and Am 2 olar ~ m 2 . — m\ = 
(0.9-1.5) x 10" 4 eV 2 . 
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FIG. 5: Fits to the atmospheric and solar neutrino data (procedure A) for the mixing matrix denoted 
"(2, 1) - +" (see Sec. 0). Contours are given at x 2 = 5, 10, 15, 20, 25. Sum also at 30, 35, 40, 45, 50. Left 
panels: CP conservation; right panels: CP non-conservation. 
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TABLE III: Best fits to the atmospheric and solar neutrino data from the fitting procedure A. The number 
of degrees of freedom is 18, masses are in eV. 



Solution 


mi 


m 2 


TO 3 


Xatm 


Xsolar 


X 2 




(1,1)- + 


0.0049 


0.013 


0.038 


19.2 


3.8 


23.0 


CP cons. 


same 


0.0036 


0.012 


0.037 


18.0 


3.7 


21.7 


CP non-cons. 


(1,1)-- 


0.0027 


0.011 


0.038 


17.1 


3.7 


20.8 


CP cons. 


same 


0.0036 


0.012 


0.038 


17.7 


3.7 


21.4 


CP non-cons. 


(1,2)+ + 


0.0047 


0.008 


0.052 


14.9 


5.8 


20.7 


CP cons. 


same 


0.0030 


0.007 


0.052 


15.3 


5.7 


21.0 


CP non-cons. 


(1,2) + - 


0.0022 


0.006 


0.052 


16.5 


6.1 


22.6 


CP cons. 


same 


0.0030 


0.007 


0.052 


16.3 


5.9 


22.1 


CP non-cons. 


(2,1)+ + 


0.0045 


0.008 


0.052 


15.1 


6.1 


21.2 


CP cons. 


same 


0.0028 


0.006 


0.052 


15.4 


5.3 


20.8 


CP non-cons. 


(2, 1) + - 


0.0021 


0.006 


0.052 


16.7 


6.0 


22.7 


CP cons. 


same 


0.0028 


0.007 


0.052 


16.6 


5.3 


21.9 


CP non-cons. 


same 


0.0029 


0.007 


0.095 


16.3 


4.5 


20.8 


CP non-cons. 


(2,1)- + 


0.0045 


0.013 


0.052 


14.5 


3.6 


18.1 


CP cons. 


same 


0.0035 


0.011 


0.053 


13.8 


3.6 


17.3 


CP non-cons. 


(2,1)-- 


0.0026 


0.010 


0.052 


12.5 


3.8 


16.3 


CP cons. 


same 


0.0034 


0.011 


0.052 


13.1 


3.7 


16.8 


CP non-cons. 



TABLE IV: Best fits to the atmospheric and solar neutrino data from the fitting procedure B. The number 
of degrees of freedom is 16, masses are in eV. 



Solution 


mi 


m 2 


m 3 


Xatm 


Xsolar 


x 2 




(2,1)- + 


0.0047 


0.011 


0.054 


16.2 


0.8 


17.0 


CP cons. 


same 


0.0034 


0.010 


0.054 


15.0 


0.8 


15.8 


CP non-cons. 


(2, 1) - 


0.0025 


0.010 


0.054 


15.0 


0.8 


15.8 


CP cons. 


same 


0.0035 


0.010 


0.054 


15.3 


0.8 


16.1 


CP non-cons. 



B. Impact of CHOOZ Data 



The CHOOZ data [ |3"4"| ] are known to disfavor models with "large" values of \U e3 \. In our 
model, large \U e3 \ typically require m 2 comparable with m 3 , as can be seen in the example shown 
in Fig. |j. Thus, we do not expect the minima shown in Table |l| to be significantly altered by 
the inclusion of the CHOOZ data. This is in fact the case, as shown in Table [V|. For the fits of 
Table [En], the inclusion of the CHOOZ data (14 data points) increases the total x 2 by 3-4 units, 
without changing the best-fit mass values. 

We show in Table |V| the relevant mixing elements corresponding to these best fits. Since 
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TABLE V: Best fits from the fitting procedure A, including the CHOOZ data [34]. Masses are in eV. 



Solution 


mi 


m 2 


m 3 


Xatm 


-^solar 


XCHOOZ 


x 2 




(2,1)- + 


0.0044 


0.013 


0.052 


14.5 


3.6 


2.9 


21.0 


CP cons. 


same 


0.0034 


0.011 


0.052 


13.6 


3.7 


3.1 


20.4 


CP non-cons. 


(2,1)- - 


0.0026 


0.010 


0.053 


12.6 


3.8 


3.9 


20.2 


CP cons. 


same 


0.0034 


0.011 


0.052 


13.1 


3.6 


3.1 


19.9 


CP non-cons. 



TABLE VI: Mixing matrix elements U a j (a = e, /i, r) for the best fits, given as (modulus, phase/7r). Masses 
are in eV. 



Solution 


mi 


m 2 


m 3 


U a l 


u a2 


u a3 


Jcp 










(0.87,0.00) 


(0.48, 1.00) 


(0.13,0.00) 




(2,1)- + 


0.0044 


0.013 


0.052 


(0.39,0.00) 
(0.31,1.00) 


(0.81,0.00) 
(0.33,1.00) 


(0.44,0.00) 
(0.89,0.00) 


0.00 










(0.86, -0.49) 


(0.50,0.47) 


(0.08,-0.38) 




(2,1)- + 


0.0034 


0.011 


0.052 


(0.43,0.04) 
(0.28, 1.00) 


(0.79,-0.01) 
(0.34, 1.00) 


(0.44,0.00) 
(0.90,0.00) 


0.0060 










(0.85,0.00) 


(0.52,0.00) 


(0.03,1.00) 




(2,1)-- 


0.0026 


0.010 


0.053 


(0.46, 1.00) 
(0.26,0.00) 


(0.77,0.00) 
(0.36,1.00) 


(0.44,0.00) 
(0.90,0.00) 


0.00 










(0.86,-0.51) 


(0.51,-0.47) 


(0.08,0.38) 




(2,1)-- 


0.0034 


0.011 


0.052 


(0.43,0.96) 
(0.28,0.00) 


(0.79,0.01) 
(0.34, 1.00) 


(0.43,0.00) 
(0.90,0.00) 


-0.0060 



these cases all refer to Solution 1 for the charged leptons, then to a good approximation we have 
U ~ R(£). Hence, the cases presented here are pairwise related (d — > —d) by the symmetry of 
R(v), discussed in Sec. |V[ 

In the decoupling approximation (neglecting f/ e3 , which ranges from 0.03 to 0.13), the rel- 
evant quantities for solar and atmospheric neutrinos are |C/ e iC/ e2 | — 0.42-0.44 and \U^U T ^\ ~ 
0.39-0.40, respectively, corresponding to large mixing in both cases. As shown in Table |VI[ when 
CP is not conserved, the amount is given by 

J CP = ±6.0 x 10" 3 (9.1) 

for neutrinos. This is much larger than the corresponding quantity for quarks as given by Eq. (|2~^). 



X. DISCUSSION 

Our theoretical description of neutrino oscillations suffers from the fact that the calculation 
of the three-flavor MSW effect relies on a somewhat crude solar model electron density. Nev- 
ertheless, it is most rewarding to find that the determination of the three neutrino masses can be 
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carried out successfully, i.e., in very reasonable agreement with the experimental data. We have 
found excellent fits to the data with m 3 = (52-54) x 1(T 3 eV, m 2 = (10-13) x 1CT 3 eV, and 
mi = (2-5) x lO" 3 eV. 
In a regime where 

mi < m 2 <C m 3 , (10.1) 

the atmospheric neutrino data determine m 3 = 0{-\/ Am^,) ~ 0.05 eV, and the solar neutrino 
data determine m 2 = (9(>/Am 2 olar ) ~ 0.01 eV. The test of the model then lies in (i) reproducing 
the hierarchy Q10.1D and (ii) the determination of m 1 . 

The fact that large mixing is required by the atmospheric as well as the solar neutrino data 
essentially forces the model into a region of parameter space where there is a strong hierarchy. To 
some extent, this can be read off from Figs. 

Another issue is to what extent the data can distinguish the different discrete parameters of the 
model, Solution 1 vs. Solution 2, as well as the b and d "parities." The data favor Solution 2 for 
the neutrino mixing and Solution 1 for the charged leptons, both by a margin of 4.4 units of x 2 . 

The best fits have been found for charged lepton Solution 1. Since the charged lepton masses 
are strongly hierarchical, the rotation matrix corresponding to Solution 1 is very close to the unit 
matrix. Thus, the overall neutrino mixing matrix U is rather close to R(v). This explains why the 
mass values obtained are close to those presented earlier JT8|]. 

It is of some interest to compare in more detail with the case 

U = R{v). (10.2) 

For the two best fits, Solutions (2) — h and (2) , we find \ 2 = 17.2 and 16.7, respectively, 

with m 3 = 0.057 and 0.052 eV. The former solution is the best fit of [Jl~8|] . 

Actually, in the limit of no mixing in the charged-lepton sector, R(£) = 1, when U = R(u), we 
see that the atmospheric transition probability (|6.3[) is invariant under b — > —b, as well as under 



d — > —d or d — > id. Similarly, for the MSW equation (j7T2|), the sign change b — > —b can be 
compensated for by a wave function sign change, 3 — > — 3 , and the sign change d — > —d (or 
d — > id) can be compensated for by sign (or phase) changes of 2 and 3 . Thus, these are exact 
symmetries of |0i| 2 in this limit of R(£) = 1. They are therefore approximate symmetries for 
Solution 1 for the charged-lepton sector. 

It is also interesting to compare with the rather different model proposed by []33|]. The masses 
they find are m 3 = 0.0506 eV, m 2 = (7.46-7.48) x 10~ 3 eV, and m x = (2.39-2.43) x 10" 3 eV. As 
mentioned above, the fact that m 3 and m 2 are rather similar to the values we find is unavoidable 
within a hierarchical fit. 

We add a comment on CP non-conservation in the lepton sector. In the work of Lehmann et 
al. [[]]], given the other parameters, CP non-conservation is maximal. For this reason, when e = i, 
the lepton CP non-conservation is also maximal in the same sense. In other words, the lepton CP 
non-conservation is either zero (e = 1) or maximal (e = i). Thus there are only two cases instead 
of a continuum. 

The absolute values of the elements of the mixing matrices are somewhat different in these 
two cases. These differences are not sufficiently large to decide experimentally whether CP is 
conserved or not in the lepton sector, since both give comparable fits, as shown in Table 0. Note 
that the x 2 is slightly smaller (by 0.3-0.6) for the case of CP non-conservation. That the Jarlskog 



determinant ( pT| ) for leptons is significantly larger than that of ( |2.9D for quarks is related to the 



fact that the ratios of neutrino masses obtained here for leptons are larger than those for quarks 
(for a more general discussion, see [J3"6|l). 
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APPENDIX A 

We discuss in this appendix some elementary properties of the mass matrix of Ref. [Jl|] 




M = I —id c b | , (Al) 
\ b 

with b 2 = 8c 2 . From (P-5|), this M is diagonalized by 

M = diag(i, 1, 1) i?diag(m 1 , -m 2 ,m 3 ) i2 T diag(-i, 1, 1). (A2) 

The rotation matrix R is discussed in Sec [IV]. The case of the mass matrix without the factors of 
% and — % is entirely similar. The relations between the elements of M and the masses have been 
given in Eq. (|2~77l). 

Consider first the triangular region (|2.8D studied in Ref. []IJ]. In this region, the Si, S 2 and S 3 of 



( |277| ) satisfy 

S3 - S l S 2 = (m 2 - mi)(m 3 - m 2 )(m 3 + mi) > 0. (A3) 



Also from (|2.7[), the 33 entry of M, namely a, satisfies the cubic equation Q3.1 1| ), and the possible 



signs of a have been discussed immediately thereafter. Since c = Si — a, the corresponding cubic 
equation for c is 

9c 3 - 10S lC 2 + (S 2 + S 2 )c + (S 3 - S 1 S 2 ) = 0. (A4) 

When ( |3.1 1| ) has three real solutions, so does ( |A4[ ). Furthermore, by the inequality (|A3|), two 
of the real solutions must be positive, while the third one is negative. A comparison with what is 
known about ( B.ll| ) shows that, in the triangular region (|2T8|), 



c < for Solution 1, 

c> for Solution 2. (A5) 

When m 2 — > rri\ or m 3 — > m 2 , it follows from ( |A3| ) that S3 — = 0. In either limit, (|A3]) 
becomes a quadratic equation when c^O, 

9c 2 - 105ic + (Sf + S 2 ) = 0. (A6) 

In the limit m 2 — * m\, this has two positive roots. In the limit m 2 — ► m 3 , this has one positive and 
one negative root. Therefore, for Solution 1, c remains negative when m 2 — > mi or m 2 — ► m.3. On 
the other hand, for Solution 2, c remains positive when m 2 — > m 1 , whereas 

c -> (A7) 

when m 2 — > m 3 . 
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From the discussion after Q3.1 1| ), the interesting case is where this cubic equation has three real 
roots. This holds not only in the triangular region ( J2~8"| ) but in a larger region. The boundary of 
this larger region is given by the straight line m 1 = together with curves obtained by setting the 
discriminant of the cubic equation Q3.11|) to zero: 



(8Sf + S 2 ) 2 (S 2 l - 36S 2 ) - 2187S'! + 2380S 3 S 3 + 2754SiS , 2 S 3 = 0. 



(A8) 



The resulting region is only slightly larger than the triangle (|2.8|), and has been discussed in Sec. [III]. 
Let the region shown in Fig. [I] be called 1Z, while 1Z denotes the triangle (|27B"|); then 1Z consists 
of TZ together with IZi (the small region where m x > m 2 ) and 1Z 2 (the region where m 2 > m 3 ). 
The extension of flA5[ ) to TZ% and 1Z 2 gives simply 



c < for Solution 1 in all of 1Z 
c > for Solution 2 in 1Z and TZi 
c < for Solution 2 in 1Z 2 . 



(A9) 



The next property to be discussed is the behavior of the m's when both c and d are small 
compared with a: 

\c\ < a and \d\ < a, (A10) 

the relative magnitude of c and d being arbitrary. In general, by Eq. (|2.7[), the masses mi, — m 2 
and m 3 satisfy the cubic equation 



A 3 - (a + c)A 2 + (ac - 8c 2 - d 2 )A + ad 2 



0. 



In the limiting case ( |A10D , one of the solutions is 

A = m 3 ~ a, 



(All) 



(A12) 



while the other two solutions are both small. These two small solutions, mi and — m 2 , are deter- 
mined approximately by the quadratic equation 



— aA 2 + (ac — d 2 )\ + ad 2 



0. 



(A13) 



The solutions of (IA13I) are 



2a 



[ac -d 2 ± y/(ac - d 2 ) 2 + 4a 2 d 2 ] 



(A14) 



or 



m 2 



l[c+Vc 2 + 4d% 



|[-C+ y/(? + 4d?]. 



(A15) 



Note that these approximate solutions are applicable only to Solution 1, because Eq. ( |3.13[ ) implies 
that, for Solution 2, c cannot satisfy the inequality ( |A1(H ). Eq. QA1 5[ ) shows directly that c < for 
Solution 1 in the triangular region (|2.8D. 
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APPENDIX B 



For completeness we give here asymptotic formulas for the 3 Fi required for reconstructing the 
2-F 2 . These asymptotic formulas turn out to be quite useful and in particular are accurate for the 
regions of the fits discussed in Sec. |VIII| . For the wave functions of the electron neutrino, tpi of 
Eq. (|74l), the two relevant 3 Fi are (cf. [|2|]) 



l + z£, 1-ir), 1+iC 
1 + 



(27T) 



/2 ^ r(i + < / )r(^)r(i-z7y-iC / ) 



r(i + iC) 



3^1 



1 + i£, 1 ~ "7, 1 - «C 
1-zC' 

ivr/2„ ( \ 1 +^ 



r(i-<C0 



r(i + ze)r(i-^)r(z77-<0 



where 



(Bl) 



(B2) 



/*00 /» 

1= / di / rfs(l + s + st) _1 exp{-i[ty -^\nt + rjlns- (C — CO ln(l + s) 
Jo Jv 

+ Cln(l + s + st)]}, (B3) 

pea poo 

1=1 dt dsil + s + st)' 1 exp{-i{ty-£\nt + fjlns + (C-C)M 1 + s ) 
Jo Jo 

-C\n(l + s + st)]}, (B4) 



and 



Note that all these quantities in Eq. ( |B5| ) are positive because 

/ii < uj 2 < /i 2 < ^3 < A«3- 



(B5) 
(B6) 

(B7) 



The double integrals for / and / of Eqs. ( P3| ) and ([B3]) can be carried out approximately using 
the method of stationary phase, and the resulting asymptotic formulas are 



27r a. 



-1/2 



e(C-C')(*o 



1-tn) 2 C(C-C') 



+ 



(i+to)ec 

'o 



c 2 



K-l)^ (a -l) 2 
x exp{-i[t y - £lnt - (rj + (') lna - (( - £') ln(<x - 1 
+ Cln(«7 -l-*o)]} , 



cr - 1 — *o 



-1/2 



(B8) 
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and 



2tt a, 



-1/2 



^C-O^o + i + to) 2 C(C-C') , (i+*oXC 



+ 



c 2 



x exp{-i[t y - £ln* -{V- (') In (T + (C - CO ln(^o + 1) 
-Cln(a + l+to)]} , 



-1/2 



where ct and t are determined by 



c 



?7 + 



while (T and i Q are determined by 



---- 

^0 (To — 1 — ^0 

C-C' (i + <b)C 



(T - 1 (T - 1 - t 



C 



o. 
o. 



rj + 



C 

y — = — - 

to (T + l+t 

C-C (i + *o)C 



o. 
o. 



(T + 1 (T + 1 + ^0 

Note that the a here is equal to —1/ s , where s is defined in the Appendix of [|2"3|]. 



(B9) 



(BIO) 



(Bll) 



APPENDIX C 



We shall provide here some further details, beyond what was given in on the reconstruc- 
tion of and ^^ 2 ' ) from the 3-Fi's. 

We introduce subscripts i and k to label the cases (5.3)— (5.5) and (5. 12)— (5. 13), respectively, 
of Then, Eqs. (5. 12)-(5. 13) of [E§ may be rewritten as 



where 



and 



ff\z) = G lk [Y m ff\z) + Y lk2 jf(z) + Y ih3 /f (z) 

7iT(l - a 2 + a±) 



G 



Y 



r(-/3i + «i)r(-/3 2 + ai)r(-^ 3 + an) 
r(-ft + ft)r(-ft + A) 1 



T(l - oti + /3i)r(l - a 2 + /3i) sin7r(a 1 - /3x) 



etc., with «j and /?_■ given by Eqs. (5.2)-(5.5) of [E4j]. Furthermore 



{^2,^2^} — {Cji,^lj} 



Note that Yj 2 j = Yjij apart from the change of argument of sin^a^ — (3j)]. 
Define now 

fW(z) = e -^ /2 P\z), f {k \z) = e -^ /2)ak ] {k) {z). 



(CI) 
(C2) 
(C3) 

(C4) 
(C5) 
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Then, Eq. ( |CI{ ) can be rewritten as 

f?\z) = B ik jP{z) + B m ff\z) + B lk3 ff\z) 

with 



B; 



(C6) 
(C7) 



For each i, we have the following two equations: 

B m fl l \z) + B a2 ff\z) = [ff\z) - B il3 ff\z) 
B m f}%) + B m ff\z) = \ff\z) - B l23 ff \z] 



(C8) 



Here, fi (z) and fj; (z) are given in terms of 3-Fi's, whereas the fi 3 \z) are given in terms of 
2 F 2 s. These are then solved for f^\z) and f- 2 \z), from which the 2 F 2 of Eq. (5.7) are obtained. 
To obtain the physical neutrino wave functions, one has to rotate back to the fa of Eq. (J7TT]). 
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